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ABSTRACT 

The  objective  of  this  project  was  to  study  the  effect 
of  small  artificial  cavities  on  natural  convection  from  a 
horizontal  surface.   Tests  were  run  with  Feron  113.   Data 
for  heat  flux  as  a  function  of  bulk  temperature  difference 
were  carefully  obtained.   These  data  yielded  Nusselt  number 
as  a  function  of  Grashof  number  or  Rayleigh  number.   All  of 
these  were  then  compared  with  the  data  obtained  by  Hiep 
[Ref.  2].   Experimental  results  are  presented  for  the  heat 
transfer  from  horizontal  circular  disks,  with  and  without 
artificial  cavities. 

The  artificial  cavities  were  found  to  affect  the 
natural  convection  heat  transfer  from  a  horizontal  surface: 
the  heat  transfer  coefficients  increased  with  the  number  of 
cavities. 
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I.   INTRODUCTION 

A.    BACKGROUND 

It  has  been  shown  by  Duncan  [Ref.  1]  that  artificial 
cavities  can  have  a  significant  effect  on  boiling  heat 
transfer.   If  the  cavities  are  of  the  correct  size  it  is 
possible  to  induce  nucleation  at  lov/er  superheats  than  would 
otherwise  be  necessary.   Recognizing  the  enhanced  heat  trans- 
fer which  artificial  cavities  can  induce  in  the  boiling  regime 
one  is  led  to  inquire  as  to  the  effect  that  these  cavities 
might  have  on  the  natural  convection  heat  transfer  which 
precedes  boiling.   Hiep  [Ref.  2]  undertook  an  experimental 
investigation  of  the  effect  of  artificial  cavities  on  the 
natural  convection  heat  transfer  from  a  horizontal  circular 
disc.   Tests  were  run  with  water  and  Freon  113.   For  water 
O'Connor  [Ref.  3]  ran  tests  with  a  mirror  finished  test 
surface  without  the  artificial  cavities  and  plotted  the  data 
in  terms  of  dimensi onl ess  parameters,  Nusselt  number  versus 
Grashof  number.   Hiep  started  with  seven  artificial  cavities 
for  water,  and  gave  the  correlation  result  of  Nusse  1 1 number 
in  term  of  Rayleigh  number  as: 

Nu  =  0.38  (Ra)°-^^  (1) 

For  Freon  113  Hiep  ran  the  tests  with  a  blank  and  seven- 
hole  test  sections,  again  he  gave  the  results  in  terms  of 
Nusselt  number  versus  Grashof  number  and  Nusselt  number  against 
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Rayleigh  number  the  correlation  results  for  blank  test  section 
was 


u  =  0.48  (Ra) 


0.30 


(2) 


The  experimental  apparatus  that  O'Connor  [Ref.  3]  and 
Hiep  [Ref.  2]  used  was  almost  the  same.   Test  surface  had  a 
diameter  of  0.895  inches,  installed  in  the  center  of  a 
cylindrical  tank  9.00  inches  in  diameter  and  4.5  inches  high 
and  set  at  the  same  level  as  the  tank  floor  (see  Figure  3). 
A  controllable  resistance  heater  as  a  heat  source  was  housed 
in  a  copper  block  attached  to  the  bottom  of  test  section. 
A  calibrated  resistor  in  series  with  the  heater  was  used  to 
determine  accurately  the  input  power  to  the  heater.   Eight 
thermocouples,  wired  in  series  were  used  to  measure  the 
average  bulk  temperature  of  the  fluid  and  one  removable 
thermocouple  was  used  to  provide  a  check  on  this  temperature. 
To  record  the  temperature  distribution  in  the  stainless  steel 
cylinder,  four  stainless  steel  sheathed  thermocouples  were 
equally  spaced  along  the  axis  of  the  test  section.   The 
temperatures  at  these  locations  were  obtained  using  a  Dana 
digital  voltmeter. 

Since  all  the  results  and  discussion  of  Hiep's  work 
depend  on  the  accuracy  of  determination  of  the  test  surface 
and  fluid  temperature  the  calibration  of  all  thermocouples  is 
extremely  important.   Hiep  used  the  standard  thermocouple 
tables  instead  of  the  average  of  his  calibration  value.   The 
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uncertainty  value  in  Hiep's  thermocouple  reading  was  0.5°F, 
which  results  in  a  large  uncertainty  for  heat  transfer 
coefficient.   The  change  in  Nusselt  number  between  the  blank 
test  section  and  the  seven  hole  test  section  as  determined  by 
Hiep  is  within  the  experimental  uncertainty  assigned  to  the 
Nussel t  number. 

The  uncertainty  in  Hiep's  experimental  data,  therefore, 
tends  to  make  it  difficult  to  draw  a  firm  conclusion  as  to 
the  effect  which  artificial  cavities  have  on  the  natural 
convection  heat  transfer  process.   His  experiments  were 
limited  to  two  test  surfaces,  blank  and  seven  cavities.   He 
also  limited  the  range  of  surface  temperature  to  those  below 
the  saturation  temperature  of  Freon  113.   In  other  words  the 
experiment  did  not  go  into  the  boiling  phase  that  Duncan 
[Ref.  1]  obtained. 

B.    THESIS  OBJECTIVES 

The  previous  discussion  has  led  to  the  objection  of  this 
experimental  works,  that  is: 

1.  To  investigate  the  effects  of  the  drilled  cavities 
on  natural  convection  heat  transfer  in  Freon  113 
with  greater  number  of  test  surfaces  than  Hiep 
obtained. 

2.  To  obtain  results  for  the  boiling  phase  of  Freon  113 
to  see  if  it  agrees  with  the  Duncan  results. 

3.  Try  to  keep  the  uncertainty  value  of  heat  transfer 
coefficient  and  other  di mens i onl ess  parameters  as 
low  as  possible. 
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1 1 .       EXPERIMENTAL    APPARATUS 

A.  DESIGN  CONSIDERATION 

1.  The  experimental  apparatus  used  in  this  study  is 
the  same  as  had  been  used  by  Hiep  [Ref.  1].   A  complete 
description  can  be  found  in  his  thesis  and  the  thesis  of 
O'Conner  [Ref.  3].   A  schematic  diagram  of  the  experimental 
set-up  is  shov/n  in  Figure  3.   For  convenience  the  ice-bath 
used  by  Hiep  was  replaced  by  an  automatic  electronic  ice-point 
for  the  reference  thermocouple.   This  automatic  ice  junction 
was  set  to  32°F  as  reference  temperature.   A  Dana  digital 
voltmeter  that  had  been  used  to  read  and  record  the  thermo- 
couples was  replaced  by  a  Hewlett  Packard  digital  voltmeter. 

2.  A  Freon  condenser  was  also  added  to  recover  the 
Freon  that  evaporates  during  the  experiment.   This  condenser 
was  connected  to  the  top  plate  of  the  tank  as  shown  in  Figure 
2,  and  was  cooled  by  tap  water. 

B.  TEST  SECTION-HEATER  ASSEMBLY 

Three  test  sections  were  used  in  this  experimental  work: 
one  blank,  one  with  seven  holes  which  had  been  used  by  Hiep 
and  one  with  thirteen  holes  (cylindrical  drilled  cavities). 
All  three  test  surfaces  had  the  same  diameter,  that  is  .896 
inch.   Four  0.040  inch  diameter  holes  were  drilled  along  the 
length  of  each  test  piece  to  a  depth  of  0.375  inches,  separated 
axially  by  about  0.5  inch  and  90°  apart  from  each  other. 
These  holes  were  for  sheathed  thermocouples  installation  in 
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Figure  1.   Photograph  of  Complete  Apparatus 
Including  Instrumentation. 
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Figure  2.   Photograph  of  Assembled  Experimental 
Unit. 
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Figure  3.   Cross  Sectional  View  of  Tank,  Test  Section  and 
Heater  Assembly. 
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Figure  4.   Photograph  of  Blank  Test  Section 
After  the  Test. 
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Figure  5.   Photograph  of  Thi rteen-Hole  Test 
Section  After  the  Test. 
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Figure  6.   Blank  Test  Section 
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21 


Hole 

(Depth 
Di  am. 


Sheathed 

Thermocoupl e 

(diam.  =  0.041  inch) 


Test  Surface 


i^igure    8.      Thi  rteen-Hol  e    Test    Section. 
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Figure  9.   Heater  Assembly 
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each,  test  section  as  shown  in  Figure  6-8.   The  thi rteen-hol e 
test  section  had  thirteen  cylindrical  holes  of  the  same  size 
as  the  holes  drilled  on  seven-hole  test  section,  that  is 
0.0375  inch  in  depth  and  0.015  inch  in  diameter,  drilled  as 
shown  in  Figure  8.   The  heater  section  [Fig.  9]  was  the  same 
as  used  by  Hiep  [Ref.  2]. 

C.    INSTRUMENTATION 

1 .    Temperature  Measurement 

ISA  type  T  copper-constantan  thermocouples  were  used 
to  measure  the  temperature  distribution  in  stainless  steel 
test  section  in  order  to  get  as  accurate  as  possible  in 
temperature  distribution.   These  four  thermocouples  were 
sheathed  in  stainless  steel  (sheathed  diameter  =  0.035  inch) 
and  the  temperature  of  the  surface  was  found  by  using  a  linear 
extrapolating  method  (see  Appendix  II,  Fig.  18). 

The  bulk  fluid  temperature  was  determined  by  using 
eight  copper-constantan  thermocouples  wired  in  series  (four 
pairs)  manufactured  using  a  Dynatech  welder.   These  thermo- 
couples were  wired  in  series  into  an  insulated  switching  box 
so  the  reading  gave  an  average  temperature  for  four  pairs  of 
these  thermocouples.   A  single  removable  thermocouple  was 
used  to  provide  a  check  on  this  temperature.   The  sheathed 
thermocouple  leads  were  securely  fastened  to  the  bottom 
phenolic  plate  and  then  to  a  Jones  bar  arrangement.   Because 
the  sheathed  thermocouples  were  difficult  to  manipulate  and 
subject  to  twisting,  and  thus  shorting,  giving  erroneous 
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reading,  the  more  flexible  wires  were  used  for  the  bulk 
temperature  thermocouples  were  used  as  extension  wire  into 
the  switching  box.   A  reference  thermocouple  was  inserted 
into  the  automatic  ice-point  reference.   A  Hewlett  Packard 
digital  voltmeter  accurate  to  0.001  millivolt  was  used  to 
read  the  output  of  the  thermocouple  in  millivolt. 
2 .    Power  to  Heater 

A  Lambda  regulated  power  supply  was  used  to  provide 
input  voltage  to  the  heater.   A  calibrated  resistor  was  placed 
in  series  with  the  power  supply  and  heater.   This  resistor 
was  mounted  on  a  large  heat  sink  and  was  isothermal  and 
therefore  of  constant  resistance  for  all  of  the  experimental 
runs.   With  this  arrangement,  a  voltage  reading  across  the 
resistor  and  across  the  heater  were  recorded.   Knowing  the 
voltage,  and  resistance,  the  current  through  the  circuit  could 
be  determined  by  Ki rchoff current  law: 


'h  =  'r  =  hf" 


(3) 


where 


R 


the  current  across  the  heater, 

the  current  across  the  resistor 
2.013  ohm. 

voltage  across  the  resistor. 


and  the  product  of  this  current  and  voltage  across  the  heater 
gave  the  input  power,  i.e.. 


or 


'   -   hh 


P  =  V^(Vj^,R)    watts. 


(4) 
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This  value  could  be  accurately  reset  for  different  runs 
The  voltage  across  the  heater  and  across  the  calibrated 
resistor  were  determined  using  the  Hewlett-Packard  digital 
voltmeter,  with  an  accuracy  of  0.001  volts.   This  power  was 
not  used  to  determined  the  heat  flux  but  only  to  set  the 
operating  point  from  run  to  run. 
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Ill  .   EXPERIMENTAL  PROCEDURE 

A.    TEST  SURFACE  PREPARATION 

It  was  very  important  that  the  test  surface  be  free  from 
any  machine  marks  and  scratches  that  could  influence  the  heat 
transfer  other  than  the  artificial  holes.   For  this  reason 
the  surface  was  mirror  finished.   The  machine  marks  were 
removed  by  using  a  Buehler  Metallurgical  Manual  grinder;  both 
0  emery  polishing  paper  and  3/0  emery  polishing  paper  were 
used.   After  that,  the  surface  was  washed  with  detergent  and 
warm  water  to  remove  the  previous  abrasive  and  then  washed 
with  alcohol.   In  order  to  get  a  mirror  finish  the  test 
surface  was  then  wet  polished  on  a  metallurgical  wheel  covered 
with  velvet.   The  wheel  was  impregnated  with  one  micron  diamond 
dust.   The  test  surface  was  placed  face  down  on  the  wheel  and 
a  low  speed  of  rotation  was  used.   During  this  work  alcohol 
was  used  to  wet  the  surface  at  all  times.   After  completion 
of  polishing,  the  test  surface  was  washed  with  alcohol  again 
and  was  dried  with  a  dryer,  and  then  wrapped  with  special 
paper  to  prevent  contamination  and  touching.   Before  installa- 
tion the  seven-hole  and  thirteen-hole  test  sections  were  put 
in  an  ultrasonic  cleaner  for  ten  minutes,  were  washed  with 
distilled  water,  rinsed  with  acetone  and  dried. 

The  seven-hole  and  thirteen-hole  test  surfaces  were 
drilled  before  the  mirror  finish,  by  using  a  Sphinx  spirec 
pivot  drill,  0.015  inch  diameter.   The  cavities  depth  were 
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controlled  by  using  a  feeler  gauge  on  the  lathe.   A  ratio  of 
approximately  2.5  was  desired  for  cavity  depth  to  diameter, 
so  the  depth  was  about  0.040  inch;  the  location  of  seven 
holes  and  thirteen  holes  were  as  shown  in  Figures  7  and  8 
respectively. 

A  height  gauge  which  had  a  scale  readable  up  to  0.0001 
inch  was  used  to  determine  the  location  of  the  center  of  the 
thermocouple  holes  relative  to  the  test  surface.   This  gauge 
was  also  used  to  measure  the  diameter  of  the  test  surface. 

B.    APPARATUS  ASSEMBLY 

The  main  parts  of  the  apparatus  were  the  test  section, 
heater  assembly,  and  the  tank  assembly.   The  test  section- 
heater  assembly  was  securely  bolted  to  the  phenolic  base 
plate,  as  described  by  Heip  [Ref.  2]  and  once  the  insulation 
was  in  place  and  the  can  installed,  the  phenolic  plate  that 
served  as  the  bottom  of  the  tank  was  carefully  fitted  to  the 
test  section.   The  level  of  the  test  section  was  adjusted 
from  the  bottom  base  of  the  assembly  by  screws.   The  phenolic 
bottom  plate  of  the  tank  was  secured  by  three  threaded  rods 
120°  apart,  with  the  phenolic  base  plate.   These  three  rods 
were  also  used  to  adjust  the  test  surface  in  level  with  the 
phenolic  bottom  plate.   For  the  tank  assembly  four  "0"  rings 
of  9.0  inch  diameter  were  used  for  sealing  to  make  sure  that 
the  cylindrical  enclosure  tank  was  leak  proof.   Two  "0"  rings 
for  each  side  (top  and  bottom)  of  the  tank  were  used.   The 
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cylindrical  plexiglass  container  was  then  put  in  place  and 
the  top  installed,  six  threaded  rods,  60°  apart  were  used  to 
secure  the  top  plate  to  the  bottom  plate  of  the  tank.   The 
glass  laboratory  condenser  was  installed  on  top  of  the  test 
tank  to  recover  any  vaporized  Freon.   Two  of  the  threaded 
rods  were  18  inches  long  so  they  could  be  used  to  hold  the 
condenser  and  insure  that  it  would  sit  tightly  on  the  1  inch 
hole  at  the  top  plate  of  the  tank.   This  1  inch  hole  was  used 
to  pour  the  fluid  into  or  out  of  the  tank.   To  hold  the 
condenser  an  adjustable  holder  was  connected  from  each  of 
the  18  inch  rods  to  make  the  condenser  stand  still.   The 
condenser  was  held  tightly  on  the  top  but  was  also  removable 
to  allow  filling  of  the  tank. 

Six  small  holes  were  drilled  through  the  plexiglass  top 
of  the  tank  to  insert  four  pairs  of  thermocouples  in  the 
tank  to  measure  the  average  bulk  temperature  of  the  fluid, 
for  each  pair  one  of  the  thermocouples  was  put  near  the 
bottom  of  the  tank  and  the  other  one  was  placed  near  the  top 
so  that  the  average  of  fluid  temperature  all  over  the  tank 
could  be  measured.   Two  holes,  one  in  the  middle  of  top  plate 
and  the  other  on  the  side  of  top  plate,  were  used  to  insert 
one  removable  thermocouple  to  the  middle  of  the  tank  to 
provide  a  check  of  the  bulk  temperature  read  by  the  other 
four  pairs. 
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C.  "   TEST  PROCEDURE 

1  .    Temperature 

A  complete  data  run  consisted  of  starting  at  the 
lowest  heater  setting,  waiting  until  steady  state  was  reached 
(usually  about  ten  hours  for  each  data  point).   The  data  was 
recorded  every    hour  during  the  daytime,  and  at  longer  inter- 
vals during  the  night.   After  steady  state  was  reached  and 
the  data  for  that  power  setting  was  recorded,  the  power  to 
the  heater  was  increased  by  approximately  5  watts.   This 
procedure  was  repeated  until  boiling  occurred  from  the  surface 
The  power  was  then  decreased  in  5  watt  increments  and  the  data 
was  recorded  after  steady  state  was  reached.   This  procedure 
was  repeated  until  the  starting  point  was  reached.   The 
procedure  of  taking  data  during  both  a  heating  and  cooling 
process  was  done  to  check  on  any  hysteresis  effects  which 
might  exist. 

The  data  recorded  for  each  setting  were  the  bulk 
temperature,  the  temperature  of  the  individual  sheathed 
thermocouples  in  the  test  section,  and  the  voltages  across 
heater  and  calibrated  resistor.   When  the  bulk  temperature 
and  the  temperature  of  the  test  section  changed  by  less  than 
half  degree  F  per  hour,  steady  state  was  assumed  to  be 
reached . 
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2 .    Determination  of  Nusselt,  Grashof  and  Rayleigh 
Numbers 

For  each  run,  at  all  input  settings,  the  temperature 
distribution  in  the  stainless  steel  test  section  are  plotted 
and  found  to  be  linear.   The  surface  temperature  could  be 
determined  within  0.5°F  from  the  plot  of  temperature  versus 
thermocouple  location  by  extrapolation,  as  shown  in  Appendix 
B.   The  properties  of  the  fluid  were  determined  at  the  film 
temperature  which  is  the  average  of  the  bulk  temperature 
(T.)  and  the  test  surface  temperature  (Tc~)  i.e., 


T 


^b  '   T3 


The  thermal  conductivity  of  304  stainless  steel 
(K33)  was  calculated  using  the  data  from  Touloukian  [Ref.  10] 
at  the  average  of  the  temperature  of  the  two  nearest  thermo- 
couples to  the  test  surface  i.e.. 


'2   '   h 


SS 


This  corresponds  to  the  average  temperature  of  the 
length  6x.   A  simple  energy  balance  at  the  surface  using 
Fourier's  Law  of  Cooling  gives: 
6T0 


K 


SS  -6^-    ^^h-\^ 


h6T 


SS 


6X 


61. 


TsTTI 


!TU/hr-ft^-°F 


(5) 
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The  Nusselt  number  was  obtained  from  the  definition 


^f 


(6) 


Substituting  values  of  h  from  (3)  into  (4)  gives 
Kc-c:  •  <5T^  •  D 


The  Grashof  number  was  calculated  using  the  diameter 
of  the  test  surface  as  the  characteristic  length: 


Gr  = 


OBD^  (Ts-T,) 


(8) 


The  Rayleigh  number  is  the  product  of  the  Grashof 
number  and  Prandlt  number: 


Ra  =  Gr-Pr 


(9) 


where 


p   ^   a  _  vC2_ 


(10) 
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IV.   RESULTS  AND  DISCUSSION 

It  was  found  that  a  minimum  of  ten  hours  was  necessary 
for  any  input  setting  to  reach  steady  state.   Once  steady 
state  was  reached,  very  good  agreement  between  Nusselt  number, 
Grashof  number  and  Rayleigh  number,  based  on  test  surface 
diameter,  obtained  during  the  heating  and  cooling  phase  of 
the  run  resulted.   The  uncertainty  bands  are  shown  on  the  plot 
of  each  curve.   The  investigation  for  the  blank  test  section 

covers  a  range  of  Rayleigh  number  from  30x10   to  350x10  ,  for 

f.  f.  p. 

seven-hole  test  section  from  50x10   to  250x10   and  25x10   to 

350x10   for  thi rteen-hol e  test  section.   A  summary  of  the 

results  for  various  data  runs  are  shown  in  Tables  I,  II,  and 

III.   For  each  data  run  the  Prandtl  number,  Nusselt  number, 

Grashof  number  and  Rayleigh  number  were  calculated.   These 

dimensi onl ess  parameters  were  determined  using  the  test 

section  surface  diameter  as  the  characteristic  length.   As 

described  in  the  experimental  procedure,  data  were  taken  at 

each  setting  as  the  apparatus  was  heated  with  increasing  input 

and  then  again  as  the  inputs  were  decreased,  cooling  the  test 

surface.   The  cooling  phase  data  are  denoted  in  the  table. 

Surface  temperature  of  the  test  section  was  determined 

by  linear  extrapolation  from  all  four  sheathed  thermocouple 

reading.   This  gave  the  surface  temperature  with  an  error  of 

approximately  . 5°F . 
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A  log-log  plot  of  Nusselt  number  versus  Grashof  number 
and  Rayleigh  number  based  on  test  surface  diameter  is  shown 
for  each  test  surface  in  Figures  10  to  17.  The  curves  show 
that  the  results  fall  into  two  regimes: 

1.  Natural  convection  which  is  the  purpose  of  this 
investigation.   The  Rayleigh  number  for  this  range  falls  in 
between  30x10   to  180x10   for  these  three  test  surface  results 
The  curves  show  the  straight  line  and  almost  have  the  same 
slope  (see  Figures  16  and  17).   The  correlation  which  results 
for  the  natural  convection  range  of  each  test  surface  with 
Freon  113  as  the  fluid  are: 

Nu  =  0.48(Ra)°*^°  30x1 0^<Ra<l 80x1 0^  for  blank  test 

section. 

Nu  =  0.53(Ra)°--^°  30x1 0^<Ra<l  80x1 0^  for  seven-hole 

test  section. 

Nu  =  0.57(Ra)°-'^°  30x1 0^<Ra<l  80x1 0^  for  thirteen- 

hole  test  section. 

2.  Boiling  phase,  the  results  of  which  must  be 
considered  with  some  caution.   This  is  because  the  bulk 
temperature  of  the  fluid  was  always  below  saturation  and 
therefore  the  boiling  washighly  sub-cooled.   In  addition  when 
nucleation  occurred  it  began  at  the  edge  of  the  test  surface 
in  the  space  between  the  test  section  and  the  phenolic  plate. 
Subsequently  nucleation  did  take  place  on  the  test  surface 

as  can  be  seen  by  the  cavitation  damage  visible  in  Figures 
4  and  5.   The  Rayleigh  number  for  this  range  started  from 
200x10   and  up.   The  curve  in  this  boiling  phase  is  shown  as 
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a.straight  line  with  slope  greater  than  the  natural  convection 
line  (also  see  Figure  16  and  17). 

In  the  natural  convection  range  the  results  show  that  the 
heat  transfer  coefficient  of  the  fluid  increases.   This  is 
consistant  with  Duncan's  result  but  disagrees  with  Hiep's. 
Hiep  stated  that  as  the  number  of  drilled  cavities  increase 
the  heat  transfer  coefficient  decrease.   The  discrepancy  with 
Hiep's  result  can  be  attributed  to  several  factors: 

1.  The  uncertainty  in  thermocouple  calibration  which 
Hiep  used  was  0.5°F  which  led  to  a  large  uncertainty  value  in 
surrace  temperature.   Surface  temperature  was  determined  by 
linear  extrapolation  from  two  sheathed  thermocouples  near  the 
surface.   This  has  a  large  effect  on  the  results  which  depend 
on  the  surface  temperature.   For  the  present  experiment  the 
value  of  the  uncertainty  in  thermocouple  calibration  was 
reduced  by  carefully  calibrating  the  thermocouples  and  using 
the  average  value  of  the  calibration.   This  reduced  the 
uncertainty  in  the  thermocouple  measurement  to  0.2°F,  which 
results  in  a  smaller  error  in  heat  transfer  coefficient. 

2.  The  number  of  test  surfaces  that  Hiep  used  (blank 

and  seven-hole)  was  insufficient  to  conclude  that  heat 

transfer  coefficient  decreases  while  the  number  of  cavities 

increases.   That  is  Hiep's  results  for  the  seven-hole  test 

section  fell  within  the  uncertainty  band  of  the  blank  test 

section.   The  present  experimental  work  used  a  thi rteen-hol e 

test  section  in  addition  to  the  seven-hole  test  section. 

These  results  show  a  consistent  trend  of  increasing  Nusselt 
number  with  increasing  number  of  cavities. 
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V.   CONCLUSIOI 


The  results  of  this  experimental  work  lead  to  the 
following  conclusions: 

1.  The  artificial  cylindrical  drilled  cavities  affect 
the  natural  convection  by  increasing  the  heat  transfer 
coefficient  of  fluid  as  shown  in  Figure  16  for  the  log-log 
plot  of  Nusselt  number  versus  Grashof  number  and  Figure  17 
for  the  log-log  plot  of  Nusselt  number  versus  Rayleigh  number. 

2.  The  correlation  results  for  this  experimental  work 
(the  horizontal  circular  heated  disk  in  cylindrical  enclosure) 
are  : 

Nu  =  0.48(Ra)"^°,  for  blank  test  section. 

30 
Nu  =  0.53(Ra)*   ,  for  seven-hole  test  section. 

30 
Nu  =  0.57(Ra)'   ,  for  thirteen-hole  test  section. 

All  of  these  correlations  were  in  the  range  of  Rayleigh 
number  from  30x10   to  180x10  .   For  the  value  of  Rayleigh 
number  from  200x10   to  350x10   the  slope  is  increased  as  seen 
in  Fi  gures  10-17. 

3.  The  log-log  plot  of  Nusselt  number  versus  Grashof 
number  and  Nusselt  number  versus  Rayleigh  number  shows  that 
the  curves  were  straight  lines  and  had  nearly  the  same  slope 
in  the  range  of  Grashof  number  from  3x10   to  20x10   and 

r  r 

Rayleigh  number  from  30x10   to  180x10  .   Between  the  range  of 
Rayleigh  number  from  180x10   to  200x10   the  slope  started  to 


47 


change  for  each  curve,  and  then  the  curves  were  straight  lines 
again  for  Rayleigh  number  greater  than  200x10   as  shown  in 
Figuresl6andl7. 

4.  Determining  the  surface  temperature  by  the  method 
of  linear  extrapolation  as  shown  in  Figure  18,  gives  the 
results  within  an  error  of  0.5°F.  This  is  insufficient 
because  the  results  show  that  if  the  surface  temperature  is 
changed  by  only  .2-.3°F,  the  value  of  Nusselt  number  varies 
by  about  +  10/^  in  the  lower  power  input  setting.  This  made 
the  log  of  Nusselt  number  change  by  about  +_12%. 
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APPENDIX  A 
THERMOCOUPLE  CALIBRATION  PROCEDURE 

The  accurate  determination  of  the  test  section  and  fluid 
temperatures  was  a  major  part  of  this  experiment  since  all  the 
results  and  discussion  are  based  on  these  temperatures.   For 
this  reason,  precise  calibration  of  the  sheathed  and  bulk 
temperature  thermocouples  was  necessary.   When  the  calibration 
was  accomplished  all  the  wiring  and  recording  instruments  that 
would  be  used  in  the  experiment  were  not  changed. 

A  Rosemount  Calibration  System,  with  a  constant  tempera- 
ture oil  bath,  was  used  for  the  calibration.   The  ice  junction 
thermocouple  which  was  sealed  in  a  glass  tube  filled  with  oil 
was  placed  in  the  electronic  ice  point  reference  junction. 
The  four  sheathed  thermocouples  and  nine  bulk  temperature 
thermocouples  were  suspended  several  inches  into  the  oil  bath. 
A  Platinum  Resistance  Thermometer  in  conjunction  with  a 
commutation  bridge  was  used  as  a  standard.   The  calibration 
was  conducted  over  a  range  from  80°F  to  220°F  and  a  maximum 
error  for  the  thermometer  for  this  range  of  temperature  was 
0.004  millivolt. 

The  results  of  these  thermocouples  reading  compared  with 
standard  thermocouple  tables.   The  deviation  from  the  table 
values  ranged  from  +.004  to  +.016  millivolt  for  sheathed 
thermocouple,  and  the  deviation  from  table  value  ranged  from 
+.006  to  +.016  millivolt  for  bulk  temperature  thermocouples. 
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Since  all  the  deviation  is  in  the  position  range  the  decision 
was  made  to  use  the  reading  value  and  subtract  the  average 
difference  from  the  table  value.   The  values  of  uncertainty 
for  sheathed  thermocouples  was  0.2°F  and  uncertainty  of  0.2°F 
for  the  bulk  temperature. 
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APPENDIX  B 
SAMPLE  CALCULATIONS 

For  Thi rteen-Hole  Test  Section: 

Input  Power  =  3.382  watts 
(Data  recorded  6  October  1975) 

A.  DETERMINATION  OF  THE  SURFACE  TEMPERATURE 

A  graphical  method  was  used  to  determine  the  surface 
temperature.   For  each  input  setting  a  plot  of  the  temperature 
versus  the  thermocouple  distance  from  the  surface,  as  shown 
in  Figure  18,  was  made  using  all  four  reading  of  sheathed 
thermocouples.   The  plots  were  linear  and  the  surface  tempera- 
ture could  be  obtained  within  an  error  of  0.5°F. 

B.  CALCULATION  OF  THE  PRANDLE  NUMBER 
From  the  definition: 


Pr  =  —  =  — 


C 


C 


(10) 


where  fluid  properties  were  evaluated  at  film  temperature 
[T 


^  :  -A_ D  ^  86.6°F] 


Fluid  properties  at  T.  =  86.6°F  [data  from  Table  II  from 
Thermophysical  Properties  of  Refrigerants  [Ref.  4]]. 
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T[°F] 
120 


115 


110 


105 


100 


95 
93.8 

90 


Power  Input  3.382  watts 
Bulk  TemperaturG  79.4°F 


± 


1 


,2    .4    .6    .8    1.0   1.2   1.4   1.6   1  .8   2.0 

Distance  from  Test  Surface,  Inches. 

Figure  18.   Determination  of  Surface  Temperature 
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1.519  Ib^/ft-hr 


04277  BTU/hr-ft-°F 

230  BTU/lb  -°F 
m 


Using  equation  (10),  we  have 


Pr  -  (^'519)(.230)  _  ^  ,. 
Kr       .04277    "  ^'^^ 
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C.    CALCULATION  OF  THE  NUSSELT  NUMBER 
Using  the  equations: 


and 


h  = 


TsxKsTJ 


(5) 


(6) 


The  value  of  K(,(.  was  determined  for  an  average  tempera- 
ture of  cylinder  over  the  length  6x.   The  average  of  the  two 
nearest  temperatures  to  the  test  surface  in  the  stainless 
steel  cylinder  were  used.   The  value  of  K^-^.  was  obtained  from 
a  plot  of  temperature  versus  K^^  based  on  values  from 
Touloukian  [Ref.  5].   Again  fluid  properties  were  evaluated 
at  the  film  temperature  using  the  data  from  Table  II  and 
Table  13  from  Thermophysi cal  Properties  of  Refrigerant  113 
[Ref.  4]  and  the  following  data  were  obtained  for  3.382  watts 
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K 


SS 


D 

61 
6T, 
6X 
P 


105. 9°F 

100.9°F 

8.926  BTU/hr-ft-°F 

.07467  ft. 

14.6°F 

5.00°F 

.04175  ft. 


96.9  lb  /ft 
m 


The  value  of  3  were  calculated  by  definition  from  the 
thermo-dynami c  properties  [Ref.  6] 


1  dp_| 
p  dT  '  p 


1  /rpQ,  I 


where 


PI+Pq 


T-j-Tq  =  Temperature  difference  of  the  given  value 
p-j  and  pq.   From  Table  13  is  2°F. 


For  thi  s  run ,  we  have 


3  =  8.668x10""^  Ij 


Using  equation  (5)  and  (6),  the  Nusselt  number  is  obtained 


h  = 


8.926(5. ) 


.04175(14.4) 


74.235  BTU/hr/ft^/°F 


^     74.235(0.07467)  ^ 

^"      0.04277        i^y.oy/ 
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D.    CALCULATION  OF  THE  GRASHOF  NUMBER 
Usi  ng  equati  on : 


Gr  =  (g6/  2)  d"*  6T, 


(8) 


The  Grashof  number  is  obtained 


Gr  =   32.174  x  8.668x10-    (0.07467)^(14.4) 
^96.9x3600^ 


8.817x10 


E.    CALCULATION  OF  THE  RAYLEIGH  NUMBER 
Using  equation: 
Ra  =  Gr- Pr 


(9) 


Thus 


Ra  =  (8.817x10*')  (8.1686) 
=  7.202x10^ 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

In  this  experimental  vyork  the  method  proposed  by  Kline 
and  McClintock  [Ref.  7]  was  used  to  obtain  the  uncertainties 
Several  assumptions  and  approximations  had  to  be  made  in 
order  to  develop  a  reasonable  estimate  of  the  uncertainties 
for  each  variable.   The  following  table  is  summary  of  these 
quanti  ties. 


Vari  abl e 
Cp 

D 

■^ss 


6X 

6 


TABLE  IV 

Uncertainty  of  Variables 

Basi  s  for  Value 

assumption  based  on  table 
values  [Ref.  4J 

accuracy  of  measurement 

assumption  based  on  table 
values  [Ref.  4] 

assumption  based  on  table 
values  [Ref.  5] 

based  on  scale  of  graphical 
method 

thermocouple  calibration 

thermocouple  calibration 

accuracy  of  measurement 

assumption  based  on  table 
values  [Ref.  4J 

assumption  based  on  table 
values  [Ref.  4] 

assumption  based  on  table 
values  [Ref.  4] 


Uncertainty  (U) 


5x10"^    BTU/lb^-°F 


2x10"^    ft 


IxlO"-" 

5% 

0.5°F 

0.2°F 
0.2°F 
4.2x10' 
1x10"^ 


BTU/hr-ft-°F 


^  ft. 
1/°R 


1x10"^  Ib^/ft.^ 
m 


2x10 


-3 


lb  /ft-hr 
m 
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The  method  that  proposed  by  Kline  and  McClintock 
[Ref.  7]  for  determining  the  uncertainty  values  is  the 
second- power  equation. 

Uncertainties  in  the  temperature  differences  6J-,     and 
6T2  were  calculated  from 

".T,  =  U      =  [(U   )2  .  (u   )2]* 
I      b   D        S         b 


and 


=  [(0.5)^  +  (CZ)"^]   =  .538°F 


^6T,  =  %-T3  =  t(U,^)^  .  (U,^)^]' 


2         1    ' 

[(0.2)   +  (0.2)  ]   =  .283°F 


Where  U^   =  uncertainty  in  T^. 
's  ^ 


U 


T  =  uncertainty  in  T, 
^b  ^ 


same  as  Uj   U^  =  uncertainty  in  T^  and  T^  respectively. 


Uncertainty  in  the  convective  heat  transfer  coefficient 
(h)  was  determined  by  using 

k 


U, 


Uk   2 


^T   2 


U 


U    2  ^ 
^6T.  ^ 


-    [(K^f)   ^  (Jt/'   *  (jF)   ^  (6T7'] 


(11) 


All  values  of  uncertainties  are  obtained  from  Table  IV 
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The  NussgH  number  is  defined  as 


Nu  = 


hD 


(12) 


and  used  of  the  second  power  equation  for  uncertainty  in  this 
number  becomes 


k 


U 


2     W         1    ^ 


uk,„  =  [(-T^)  ■»■  (-#)  +  {tA)  ] 


Nu 


(13) 


It  can  be  shown  that  the  values  of  (-j^-)   and  -t- —  are 
very  small  and  can  be  neglected.   Therefore, 


Nu    h 


(14) 


To  find  the  uncertainty  of  the  Grashof  number  use  this 
rel ati  on  : 


so 


Gr  =  (g3/v^)  D"^  6T,   . 


Gr 


2    ^T.  2  ^ 


ih)'    +  (2^p)   +  (2-')  +  (3^0)"  +  (t^)^]     (15) 


D 


The  values  of  U  are  found  from  Table  IV. 
And  by  definition  for  Prandtl  number  is: 

Pr  =  r-^- 
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so 


Pr 
Pr 


U   2    U   2    U   2  % 


against  all  the  values  of  U  are  found  from  Table  IV 
the  uncertainty  value  for  the  Rayleigh  number, 
where 

Ra  =  Gr-Pr 


7  ?  > 

Un         "^r  ^n 

Ra    '-^  Gr  ^   ^  ^  Pj.  ^  J 


(16) 
Find 


(17) 


A  summary  results  of  this  analysis  are  given  for  each 
run  in  Table  V  for  blank  test  section  and  in  Table  VI  for 
seven-hole  test  section  and  for  thi rteen-hol e  in  Table  VII 
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TABLE  V 
Uncertainties  in  the  Experiment  (Blank) 


Percent  Uncertainty 

Grashof        Rayleigh 
Number         Number 


8.40  8.70 

4.04  4.67 

3.10  3.90 

2.20  3.30 

2.16  3.23 

2.08  3.19 

1 .99  3.14 

1.89  3.13 

1 .83  3.10 

1.80  3.10 


Power  In 

DUt 

(Watts 

Nussel t 

Number 

1  .167 

20.09 

3.140 

9.12 

4.503 

7.33 

7.640 

5.90 

8.440 

5.80 

9.536 

5.69 

11  .330 

5.55 

24.500 

5.30 

29.655 

5.28 

34.735 

5.26 
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TABLE  VI 
Uncertainties  in  the  Experiment  (Seven-Hole) 


Power  In 
(Watts 

put 
) 

Nussel t 
Number 

2.429 

11  .30 

■  5.300 

6.80 

7.791 

5.98 

10.060 

5.66 

14.970 

5.42 

Percent  Uncertainty 

Grashof         Rayleigh 
Number  Number 


5.06 
2.90 
2.32 
2.09 
1  .95 


5.64 
3.74 
3.33 
3.19 
3.12 
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TABLE  VII 

Uncertainties  in  the  Experiment  (Thi rteen-Hol e ) 

Percent  Uncertainty 

Grashof  Rayleigh 

Number  Number 

10.26  11.09 

4.00  4.64 

2.75  3.64 

2.15  3.22 

1 .98  3.14 

1.90  3.13 

1.91  3.16 


Power  In 
(Watts 

5 

ut 

Nussel t 
Number 

1  .167 

23.20 

3.382 

8.48 

5.758 

6.42 

9.176 

5.68 

16.300 

5.39 

26.280 

5.31 

34.  85 

5.29 
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